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SIMULATIONS  OF  AXISYMMETRIC  EROSION  IN 
IFR-TRANSPORTED  BEAMS 


I .  Introduction 

An  electron  beam  propagating  along  a  low  density  plasma  channel  in  the 

ion  focused  regime  (IFR)  is  characterized  by  a  freely  expanding  head,  a 

pinch  region,  where  the  plasma  electrons  are  ejected  from  the  channel  and 

the  beam  pinches  to  its  equilibrium  value,  and  the  beam  body.  When  Nc  < 

,  where  is  the  line  density  of  the  plasma  channel  electrons  and  is 

the  line  density  of  the  beam  electrons,  the  channel  electrons  are  fully 

expelled  from  the  vicinity  of  the  beam  and  the  remaining  ions  reduce  its 

space  charge  potential.  This  has  allowed  stable  propagation  of  intense 

beams  over  distances  of  several  tens  of  meters,  as  demonstrated,  for 

example,  in  the  Advanced  Test  Accelerator . ^  The  pinch  region  of  the  beam 

is  associated  with  increasing  net  current  and  azimuthal  magnetic  field  and 

an  inductive  E^  spike.  Energy  loss,  due  to  the  inductive  field,  causes 

electrons  to  be  removed  from  the  beam  head.  Additional  particle  losses  can 

result  from  the  finite  emittance  of  the  portion  of  the  beam  that  is  not 

fully  pinched.  Both  of  these  processes  cause  the  pinch  point  to  erode 

backwards  in  the  beam  frame.  These  beams  are  also  subject  to 

instabilities,  such  as  the  transverse  two-stream  and  ion  hose 
2  3  4 

instabilities.  ’  ’  In  the  present  paper,  we  work  within  the  constraint  of 
azimuthal  symmetry,  in  which  the  dominant  (transverse)  instabili ties  are 
not  present. 

The  physics  of  TFR  piopagation  lias  developed  over  a  number  of  years 

2  5 

and  was'  recently  discussed  in  some  detail.  We  will  focus  on  the 

transport  losses  associated  with  inductive  and  emittance  effects  at  the 
beam  hear).  To  complement  and  illuminate  the  previous  analyses,  we  perform 
particle  simulations  in  sufficient  detail  that  erosion  rates  can  be 
measured  directly  from  the  motion  of  the  pinch  point  relative  to  the  motion 
of  the  beam  body.  Such  simulations  have  not  been  previously  presented. 
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The  simulations  are  a  product  of  FRIEZR,  a  fully  electromagnetic  two- 
dimensional  (r,z)  particle  code  that  has  recently  been  developed  to  study 
IFR  propagation  and  transport  of  electron  beams.  It  features  fully 
relativistic  beam  and  channel  electrons  and  treats  the  ions  as  an  immobile 
fluid.  A  notable  feature  of  the  code  is  that  it  solves  Maxwell's  equations 
and  propagates  the  beam  in  a  coordinate  system  which  moves  at  the  speed  of 
light.  The  coordinates  are  r,C,T  where  C  =  ct  -  z  and  t  =  t.  An  eroding 
beam  appears  to  move  slowly  backwards  in  these  coordinates.  FRIEZR  is 
discussed  in  the  appendix. 

The  remainder  of  this  paper  proceeds  as  follows.  In  section  II  we 
derive  a  new  expression  for  steady-state  inductive  erosion  which  is  valid 
in  cases  where  scattering  and  secondary  ionization  may  be  neglected. 
Simulation  results  are  then  presented  for  comparison.  In  section  III,  we 
cite  previous  analyses  of  emi t tance-driven  erosion  and  compare  these 
results  with  simulations.  Here,  we  will  find  that  the  assumption  of  self¬ 
similar  radial  expansion  of  the  beam  head  that  is  contained  in  these 
analyses  does  not  hold.  The  simulations  provide  some  insights  into  the 
relevant  physics.  Section  IV  concludes. 


II .  Inductive  Erosion 

In  IFR  propagation,  where  energy  losses  due  to  scattering  may  be 
neglected,  the  inductive  (ohmic)  erosion  rate  per  unit  propagation  distance 
of  a  beam  pulse  may  be  simply  approximated  by 

jo _ _£  _  AC  max  . ,  . 

gp  ~  Az  =  range  ’ 


where 

6p  is 


C  is  the  pulse  length,  (3  is 
max  r  o 

the  velocity  of  the  pinch  point 


the  beam  velocity  normalized  to  c  and 

so  that  (3  -  6  is  the  rate  of  beam 

o  Kp 


loss  with  respect  to  time.  Note  that  the  first  equality  above  implies  that 

losses  are  confined  to  the  beam  head,  as  is  the  case  when  scattering  and 

radiation  losses  may  be  neglected.  This  also  requires  that  the  beam  be 

matched  to  the  plasma  channel,  so  that  the  beam  body  is  in  equilibrium. 

The  range,  which  appears  in  Eq .  (1)  above,  may  be  calculated  from  the 

beam  energy  that  is  lost  to  the  inductive  axial  electric  field  per  unit 

propagation  distance  and  the  available  energy  of  the  beam  electrons.  Using 

a  reduced  set  of  Maxwell's  equations  for  relativistic  beams,  the  electric 

2 

field  at  a  fixed  radius  r  =  r  has  been  calculated  to  be 

c 

Er<vc>  - 1  Mr)  fr(£V  ■  <2> 


where  r  is  the  radius  of  the  plasma  channel,  f  =  N  /N.  is  the  ratio  of  the 
c  c  b 

line  densities  of  channel  ions  to  beam  electrons  and  I,  is  the  beam  current 

b 

(taken  to  be  positive),  rises  from  zero  at  (  =  0  to  reach  a  constant  value 
behind  the  pinch  point.  In  deriving  (2)  it  was  assumed  that  the  plasma 
electrons  are  ejected  from  the  channel  and  move  out  to  a  conducting 
boundary  at  a  finite  radius  r  .  Beam  energy  lost  to  the  induced  field  is 
transferred  to  the  channel  electrons  as  they  are  moved  to  large  radii  and 
the  channel  potential  is  established.  The  assumption  of  a  finite 
conducting  boundary  is  necessary  because  of  the  logarithmic  naLtire  of  this 
potential.  Averaged  over  the  entire  pulse  (2)  gives 


max 


max 


F,z(rc,C)  dC 


L  fl, 
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max 
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where  L  -  21og(ry/r  )  is  a  dimensionless  inductance.  The  range  may  now  ho 
•  1 1  cm  1  a  ted  as 


range 


7 

-  Yf)mc‘ 


(A) 


here  m  is  t he  electron  mass,  (y.  .  I  )  m: ■'  is 

1  nj 
n 

;e:  tion,  and  (y,  l)mcc  is  the  final  energy 
t  from  the  beam.  Note  that  unless  tin-  beam 
.  mu  channel  ( f  =  1 )  the  propagating  gamma  o 


the  energy  of  the  beam  upon 

of  the  particles  as  they  ai< 

is  fully  neutralized  by  the 

the  beam,  y  ,  is  less  than 
o 


Yinj  '  Yo  f  YU 

2 

life  y^mc  '  is  the  field  potential  energy: 


y 


U 


eIbLb(1  -  f) 
me3 


(5) 


(6) 


Here,  L.  -  ?log(r  / r ,  )  is  a  dimensionless  inductance  and  r.  is  the  beam 
b  e  w  b  b 

i ad  ins.  Note  that  for  the  vacuum  case  (f  =  0)  the  above  gives  the  familiar 
<  .-  press  ion  for  the  limiting  current  of  a  beam  with  Y  =  Yy  propagating  in  a 
ti  il!  tube  of  radius  r  .  Using  F.qs.  (’),  (A)  and  (b),  the  erosion  rate  may 

V 

,i  itten  in  terms  of  equilibi  i  urn  boa  i  and  channel  parameters,  (3^  and 


e  f  J .  L  B 
be  p 


me  (Y, 


Yf  > 


(7) 


The  value  of  y^  can  he  estimated  via  a  simple  argument.  In  the  r;r  e 
IV) .  (/)  holds,  the  erosion  pi  <><  oss  i  •  t  a  h  I  i  shes  a  constant  velocity 
t ho  pinch  point,  v  (3  r.  Assuming  tin  enoigy  ot  the  beam  is  lo:  t  i 


t 


the  spike  adiabatically,  a  particle  with  a  velocity  infinitesimally  less 
than  Vp  will  fall  behind  the  pinch  region,  where  the  E z  spike  is  localized. 
The  final  velocity  of  this  beam  particle,  v^,  is  therefore  approximated  by 

8  =  ef  =  a  -  =  i  -  ~  (8) 

p  Yf  2y\ 


Substituting  0^  from  (8)  into  (7),  we  obtain  a  cubic  equation  which 


determines  y 


f ' 


Yf  -  (Y0  +  Yu  +  2CY2o)y2£  -  Y^ Yf  +  (Y0  +  Yy  +  C)y\  =  0  (9) 


where 


C  = 


efl.  L 
b  c 


me 


(10) 


To  verify  the  theoretical  inductive  erosion  rates  predicted  by  Eq . 

(7),  simulations  with  FRTEZR  were  performed  over  a  range  of  values  for  the 
beam  current,  1^  =  2-10  kA,  and  the  plasma  ionization  fraction,  f  =  0.2- 
1.0,  with  yq  =  10.  The  beams  generally  showed  the  characteristic  trumpet 
shape  as  expected,  but  also  had  features  that  were  not  accounted  for  in  the 
simple  model  that  produced  Eq.  (7),  such  as  the  multiple  E^  spikes  that 
were  observed  at  high  1^  and  f  values.  In  these  cases,  oscillatory 
behavior  at  the  beam  head  was  characteristic  of  the  relativistic  electron- 
electron  two-stream  instability. 

In  the  relativistic  limit,  where  the  beam  electrons  have  longitudinal 
and  transverse  effective  masses  which  differ  significantly,  the  fastest 
growing  mode  of  the  two  steam  instability  is  characterized  by  a  wave  vector 
that  is  not  parallel  to  the  direction  of  propagation.  Here,  we  consider  a 
cold  beam  interacting  with  a  cool  plasma  of  density  n^,  where  n^  is  the 


6  2 

channel  electron  density.  Typically  ’  fhe  fastest  growing  wave  vector 

_  i  2  2  - 1 

features  k.,  =  c o  /v,  and  k.  =  k  =  r.  ,  where  w  =  4tin  e  /m  and  r.  is  taken 
1 1  c  b  1  r  b  c  c  b 

to  be  the  lowest  available  k^ .  Because  we  are  constrained  by  axisymmetry, 

kQ  =  0,  giving  the  mode  a  sausage-like  structure  in  r  and  C.  Results  of  a 

FRIEZR  run  in  which  these  oscillations  were  observed  are  shown  in  Fig.  1. 

[he  plot  of  E versus  C  on  axis  shows  the  oscillations  and  corresponding 

lot  of  R  2*  the  radius  within  which  one  half  of  the  initial  beam  curient 

■■  enclosed  at  a  given  C,  shows  the  sausage  structure  of  the  instability. 

'ore  that  in  the  figures,  the  beam  is  moving  to  the  left.  Note  also  that 

.lie  beam  has  been  assigned  a  positive  chaige  so  that  E  <  0  is 

decelerating.  As  can  be  seen  in  the  figures,  the  oscillatory  behavior  is 

confined  to  the  beam  head,  where  the  channel  electrons  have  not  yet  been 

fully  ejected.  This  means  that  over  the  length  of  the  interaction  region, 

the  channel  electron  density  is  dropping  from  its  peak  initial  value  to 

zero.  In  this  r  :.r. ,  I,  --  5  \u\,  r,  -  1  cn.,  r  -  2  cm  and  f  -  1.  With  these 

b  b  c 

parameters,  the  initial  average  channel  electron  density  is  nc  =  8.3  x  1()^(* 

3 

cm  so  that  the  oscillations  will  have  an  initial  wavelength  of  =  11.6 
cm  which  increases  as  the  channel  density  drops.  The  three  distinct  peaks 
in  the  E ?  plot  of  Fig.  1  show  increasing  separation  consistent  with  this 
scenario.  Channel  electron  and  beam  electron  densities  on  axis  are  also 
plotted  in  Fig.  1.  Note  that  the  oscillations  in  F,  and  R^^  correspond  to 

oscillations  in  n  .  Note  also  that  the  n  and  n.  plots  show  the  value  on 

c  c  b 

axis,  rather  than  the  lower  average  value  in  the  channel.  Because  the 
values  of  all  parameters  are  changing  ovet  the  course  of  each  oscillation 
wavelength,  the  instability  does  not  lend  itself  to  a  more  detailed 
analysis.  ft  is  of  interest,  however,  that  this  phenomena  is  confined  to 
the  beam  head  and,  as  we  shall  see  below, 
aigument  that  produced  F,q .  (7). 


does  not  adversely  at  feet  the 


Beam  head  erosion  is  ideally  measured  from  simulation  results  by 

observing  the  movement  of  the  pinch  point,  as  a  function  of  time,  to  larger 

C  values,  where  C  -  ct  -  z.  In  the  presence  of  the  complicated  beam  head 

dynamics  discussed  above,  however,  the  location  of  the  pinch  point  is 

sometimes  difficult  to  define.  In  all  cases,  we  looked  at  a  number  of 

features,  but  the  most  reliable  measurement  has  proven  to  be  Ig  ^  versus 

C,  where  I  .  is  defined  to  be  the  current  enclosed  inside  the  radius  r  = 

5  r^.  Examples  of  this  diagnostic  are  given  in  Figs.  2  and  3.  Figure  2 

shows  superimposed  plots  of  Ig  ^  versus  C  at  ct  =  5  meters  and  10  meters 

for  the  beam  pictured  in  Fig.  1.  The  superimposed  plots  give  a  reasonable 

measurement  of  the  erosion  rate,  despite  the  dynamics  at  the  beam  head. 

Figure  3  is  a  similar  superposition  of  I  ^  versus  C  plots  at  ct  =  30  and 

40  meters.  The  parameters  are  1^  =  2  kA,  f  =  0.2,  r^  =  1  cm  and  r  =  2  cm. 

Corresponding  plots  of  ar>d  nc  versus  C  (see  Fig.  4)  show  that 

this  beam  lacks  the  instabilities  seen  in  Fig.  1,  having  a  well  defined 

pinch  point  and  E^  spike.  In  such  cases,  these  other  features  corroborated 

the  erosion  rate  estimated  from  the  I  ,  plots. 

enc  1  r 

Erosion  rates  for  various  values  of  I,  and  f  are  shown  in  Table  1.  In 

b 

all  cases,  =  10,  r^  =  1  cm,  r^  =  2  cm,  r^  =  33  cm,  (beam  rise  length) 
-  50  cm,  and  Cmax  =  100-400  cm.  The  erosion  rate  predicted  by  Eq .  (7)  is 
also  given.  These  rates  were  seen  to  be  insensitive  to  the  channel  radius, 
as  expected,  and  additional  runs  at  1^  =  10  kA  and  f  =  0.2  showed  correct 
scaling  with  Yn  and  with  r  .  It  is  interesting  that  the  results  agree  so 
closely  with  theory,  even  when  the  beam  dynamics  are  complicated.  It  is 
also  fine,  however,  that  as  long  as  the  energy  losses  are  confined  to  the 
beam  head,  ns  :r.  those  rases,  the  energy  balance  argument  that  produced  Eq . 

( / )  still  app 1 i es . 
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It  is  important  to  note  that  Eq .  (7)  differs  from  previous  analyses 

mainly  in  the  treatment  of  y^ ,  which  is  calculated  by  solving  Eq .  (9).  A 

2 

common  feature  of  previous  treatments  of  inductive  erosion  in  this  regime 

is  to  assume  either  y  >>  y,  or  y,  =  1.  It  is  easy  to  see  that,  for  our  1, 

o  t  t  J  h 

k.A  cases,  either  of  these  approximations  will  cause  the  erosion  rate  to 

-  under es t  i  ma ted  by  as  much  as  50%,  because  the  particle  dynamics  limit 

a-.ount  of  energy  that  can  be  extracted  from  the  beam.  Another 
'•r.soquer.ce  of  these  dynamics  is  that  a  population  of  particles  with  y  =  y. 
■ay  remain  behind  the  beam  head.  In  those  cases  where  y^  is  large,  as  in 
•  .o  I  t  -  7.  !*A  cases,  a  significant  portion  of  these  particles  may  remain  a; 
~.a  1  i  radii,  near  the  beam,  and  will  show  up  in  the  I  j  plots  as  an 
ap.pa:  vn t  current  enhancement  behind  the  beam  head  (see  Fig.  3).  Because 
rheso  pai  ricles  are  behind  the  pinch  region,  they  are  an  artifact  with  no 
effect  on  the  beam  head  dynamics. 

III.  Emit  tance  Driven  Eros  ion 

For  a  beam  propagating  in  near  equ i 1 i br i urn ,  transverse  emittance 
losses  are  confined  to  the  beam  head,  where  the  pinch  force  is  not 
efficient  for  confinement.  In  cases  where  this  is  the  dominant  erosion 
:cerhani sm,  the  loss  of  particles  in  the  radial  direction  causes  the  length 
r,ale  for  pinching  to  increase  so  that  an  increasing  portion  of  the  beam  is 
■■'!'■' ept  ;  ble.  One  expects  these  losses  to  be  significant  only  when  the 
pinch  for'  e  and  inductive  F,  spike  at  the  head  of  the  beam  are  weak,  as  in 
’no  'O'c  of  low  current,  When  consider  i  ng  the  doctea.se  in  the  erosion  into 
as  •  n«*  'irrerit  is  lowered,  emittance  losses  would  appeal  t  o  place  a  lowet 
so  .t.d  on  *  he  e  i  j  on  rate  gi  yen  by  F.q  .  (/). 

Ihi'  issue-  was  recently  addressed  in  this  regime  by  Buchanan,"  who 
:. .  o  ;<■'<  •  fr<-  c-n /elope  equations  of  l.oo  and  Cooper  ,  ^  which  i  equ  i  I  e  that  the 


X 


radial  profile  of  the  beam  be  a  function  that  is  parameterized  by  R(C,z), 
the  root  mean  square  radius.  This  radial  function  is  specified  and  R(C,z) 
is  solved  for.  In  Buchanan's  analysis,  a  Gaussian  profile  was  used.  The 
resulting  equations  could  not  be  solved  analytically,  but  the  lack  of  an 
equilibrium  solution  (characterized  by  8/9z  -  0)  was  demonstrated, 
indicating  a  nonzero  erosion  rate  for  all  time.  When  the  equations  were 
solved  numerically  for  a  wide  range  of  parameters  it  was  found  that  the 
asymptotic  erosion  rate  is  a  constant  given  by  empirical  expression 


0  =  0.18—  f— ] 

P  lrcJ 


2( 1+f ) 


f 


0.9 


(11) 


g 

This  was  addressed  independently  by  Nguyen  and  Uhm,  who  again  applied 

the  envelope  equations  of  Lee  and  Cooper  to  a  beam  with  a  Gaussian  radial 

profile.  They  made  the  additional  assumption  that,  after  a  short 

transient,  the  erosion  rate  is  constant  and  the  axial  profile  of  the  beam 

is  fixed  so  that  all  relevant  qualities  vary  only  as  X  =  v  t  -  z ,  where  v 

P  P 

is  the  velocity  of  the  eroding  pinch  point.  As  above,  the  envelope 
equations  were  solved  numerically  with  resulting  erosion  rates  summarized 
by  the  expression: 


0 


o 


0.094  f0,85 
0.044  f0,5 


f  >  0.1 
f  <  0.1 


(12) 


'•  f. .  ( 1 )  and  F.qs .  (11)  and  ( 1 ,? )  we  expect  that  omittance  erosion 
should  no*  contribute  significantly  to  losses  in  high  cut  rent  beams,  which 
is  apparently  the  case  for  the  parameters  used  in  Table  1.  At  low'  current 
(<  1  kA)  and  high  y  (>  ?0),  however,  these  equations  suggest  that 
sign i f i <  an  t  of  f  ec  t  s  should  be  observable  in  the  simulations. 


We  investigated  this  possibility  by  making  several  runs  with  1^=1 

kA,  y  =  20,  C  =  10  cm,  C  =  100-400  cm.  The  values  used  for  the 
o  r  max 

ionization  fraction,  0.055  <  f  <  1.0,  and  the  ratio  of  beam  to  channel 
radius,  1.0  <  r^/i'c  <  3.0,  were  chosen  so  as  to  satisfy  two  conditions: 

(1)  we  require  n^  <  n^  throughout  the  channel  so  that  all  channel  electrons 
viiJ  be  expelled  and  (2)  the  predicted  emittance  erosion  rate  must  dominate 
•er  the  inductive  effects.  From  Eqs.  (11)  and  (12)  erosion  rates  ranging 
■com  -  g  <  0.01  up  to  gQ  -  g^  >  0.08  were  expected  for  these 
parameters . 

As  in  the  inductive  cases,  erosion  was  measured  from  plots  of  I 

erst  s  C  once  steady-state  erosion  was  established.  In  contrast  with  the 

inductive  case,  where  this  could  be  determined  from  the  average  particle 

energy  in  the  beam  head,  we  had  to  rely  on  axial  profiles  of  beam 

quantities  to  indicate  steady-state,  or  near  steady-state  behavior.  In  all 

cases,  we  found  erosion  rates  of  .01  <  gQ  -  g^  <  .03,  which  varied  only 

sluggishly  with  f  and  rj /r  .  The  measured  rates  also  varied  with  the 

numerical  parameters,  AC  and  Ax,  indicating  that  the  erosion  rates  were 

less  than  or  of  the  order  of  the  numerical  errors. 

To  clarify  these  data,  we  repeated  several  of  the  runs  with  yq  =  5. 

This  reduction  in  yo  should  double  the  rate  of  emittance  erosion,  according 

to  (11)  and  (12).  Inductive  erosion  effects  were  removed  in  these  cases  by 

specifying  F.  -  0  in  the  part  of  the  simulation  where  the  beam  particles 

are  acted  on  by  the  fields  (the  particle  push  routine).  In  all  cases,  we 

found  g  g  -  0.005,  a  reduction  in  the  erosion  rate  that  is  opposite  the 
o  p 

scaling  predicted  in  (11)  and  (12).  Because  the  measurements  were  made  in 

the  asymptotic  limit  as  required  to  obtain  steady  state  behavior,  we 

<■  or,'  bided  'hat  the  rate  of  erosion  due  t  r>  the  finite  emittance  of  the  beam 
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tends  toward  zero  in  the  asymptotic  limit.  This  conclusion  is  supported  by 

further  analysis  of  the  simulations. 

To  further  investigate  the  results,  the  beam  was  diagnosed  for  average 

transverse  emittance  as  a  function  of  C  throughout  the  simulation.  Plots 

of  I  i  versus  C  from  a  run  in  which  y  =  5,  f  =  1,  r,  =  r  =2  cm,  I,  =  1 
end  o  b  c  ’  b 

kA  and  in  which  E  =  0  were  specified  in  the  particle  push  routine  as  shown 

in  Fig.  5.  Early  in  the  simulation,  at  ct  =  4  meters,  the  initial  rise 

length  of  10  cm  has  increased  to  100  cm,  and  a  low  emittance  population  of 

particles,  wi th  T  <  T  /2,  has  developed  at  the  beam  head,  where  T^  is  the 

initial  transverse  temperature  of  the  beam.  Beyond  ct  =  4  meters,  the 

axial  profile  evolves  slowly,  with  both  the  rise  length  and  the  low- 

emittance  population  continuing  to  increase  and  the  rate  of  increase  in 

these  quantities  tending  toward  zero.  Note  that  because  of  the  low  energy, 

y  =  5,  g  =  0.98,  the  beam  motion  relative  to  the  coordinates  will  be 
o  o 

discernible  over  the  course  of  th"-  simulation,  which  continues  until  ex  = 

40  meters.  This  motion  can  be  seen  in  the  the  ex  =  4  m  and  cx  =  40  m  plots 
of  Fig.  5.  When  measuring  the  cumulative  effects  of  the  transient  erosion 
process  over  a  long  distance,  this  must  be  taken  into  account. 

The  progress  toward  an  asymptotic  state  can  be  seen  in  Fig.  6,  where 
the  growth  of  the  low  emittance  population  is  shown  by  plotting  Cj.  versus 
z,  where  Cj.  is  the  cooidinate  of  the  point  at  which  the  beam  particles  have 
transverse  temperature  T  =  Tq/2  and  z  =  cx.  Particles  with  C  <  d  will 
have  lower  emittance.  As  the  beam  reaches  an  asymptotic  state,  the  motion 


of  Cp(z)  will  match  the  motion  of  tiie  beam  body,  t^(z)  =  (1  g)z,  so  that 
Cj.(z)  and  C^fz)  will  have  the  same  slope.  A  line  with  slope  dC/dz  =  1  -  g 


is  also  plotted  in  Fig.  6.  Note  that  for  z  >  20  m  the  two  slopes  have 
converged  to  within  numerical  fluctuations. 


The  result  of  an  asymptotically  vanishing  erosion  rate  is  clearly  a 
departure  from  the  predictions  of  Eqs.  (13)  and  (14).  The  discrepancy 
resides  'n  the  assumption  of  self-similar  radial  structure  which  is 
fundamental  to  the  derivation  of  the  envelope  equations.  This  assumption 
appears  to  be  too  limiting  in  this  case,  where  the  beam  must  be 
parameterized,  at  the  very  least,  by  T(C,z)  as  well  as  R(C>z).  The  beam 
•i:r/  in  fact  assume  a  highly  non-Gaussian  profile  at  early  times,  when  high 
femperature  particles  escape  and  low  temperature  particles  remain  trapped 
by  the  weak  pinch  potential.  These  trapped  particles  contribute  to  the 
increased  rise  length  of  the  beam. 

I V .  Conclusions 

We  have  seen  that  numerical  simulations  of  inductive  erosion  agree 
closely  with  the  theory  over  a  wide  range  of  parameters,  even  when  the  beam 
dynamics  are  complicated  by  unstable  behavior  in  the  beam  head.  These 
results  demonstrate  that  at  low  energies  (<  5  MeV)  careful  consideration  of 
the  relativistic  dynamics  of  the  beam  particles  are  required  for  accurate 
prediction  of  inductive  erosion  rates.  This,  in  effect,  limits  ones 
ability  to  reduce  the  erosion  rate  by  various  means  because  a  lower  erosion 
rate  results  in  a  higher  y^.  At  higher  energies,  the  yq  >>  y^ 

approximation  holds,  and  our  results  are  in  agreement  with  previous 

2 

f ormu  lat ions . 

In  the  case  of  emittance  erosion,  simulations  indicate  that  a  steady- 

state  erosion  rate  is  not  established  and  that  the  radial  expansion  at  the 

beam  head  is  not  self-similar,  thereby  departing  from  previous  theoretical 
2  8 

assumptions.  ’  As  the  beam  channel  equilibrium  is  established,  these 
beams  develop  a  long  rise  length  of  low  emittance  particles,  with  an 
erosion  rate  that  tends  towards  zero.  Because  the  simulations  show  that 
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this  asymptotic  behavior  establishes  itself  only  over  large  propagation 
distances,  this  may  not  be  easily  observable  in  laboratory  cases  of 
interest.  For  these  cases  a  new  analysis,  in  which  the  radial  dynamics  are 
treated  more  carefully,  would  be  required. 
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Appendix 

For  these  studies  we  have  used  an  electromagnetic  particle  simulation 

code,  FRIEZR,  which  was  written  to  study  the  propagation  of  high  energy 

electron  beams  in  the  Ion  Focus  Regime  in  an  axisymmetric  coordinate 

system.  The  elements  of  the  code  are  beam  and  plasma  electrons  which  are 

treated  as  particles,  and  plasma  ions  which  may  be  treated  as  an  immobile 

oackground .  The  code  is  written  in  a  transformed  coordinate  system 

(r,C,t),  where  C  =  ct  -  z  and  t  =  t.  For  a  highly  relativistic  beam,  the 

axial  variable,  C,  measures  distance  behind  the  beam  head.  The  beam  moves 

slightly  slower  than  c  so  the  beam  head  slips  slowly  backward  in  the 

coordinate  system,  but  for  high  energy  beams  and  short  propagation 

distances,  the  slippage  is  not  significant  relative  to  the  size  of  the 

simulation  region.  The  particle  dynamics  are  determined  from  the  solution 

of  Maxwell's  equations  written  in  the  (r,C,"t)  variables  with  current  and 

charge  densities  determined  from  the  simulation  particles.  We  write  the 

equations  for  the  vector  and  scalar  potentials,  A  and  (f>,  in  a  gauge 

9 

suggested  by  Lee  which  is  particularly  convenient  for  axisymmetry.  In 
this  gauge,  which  has  the  condition 

=  0  ’  (ai> 

we  may  write  the  parallel  and  perpendicular  components  of  Ampere's  law  as 


,  Ifl)  .  is,  t 

i  I9z  c  3tj  cl  1 


„  ,  3F, 

n2  a  1  z 

1  z  c  3t 


Vf  change  variables  to  the  (  r , C- c t -z  ,  t= t )  coordinate  system  and  write 
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(A4) 


(A5) 


(A6 ) 


( A7 ) 


(A8) 


(A9) 


(A10) 


(All) 


If) 


In  axisymmetry,  A£ 
<f>  equations,  so  that 


V2A 

e  2 

r 


=  0  and  the  Afl  equation  is  decoupled  from  the  a  and 


1  L_d  <L  91_)a  _ 
c  a-clc  3t  +  ad  e  "  c  e  • 


(A12 ) 


These  equations  determine  the  dynamics  of  the  system  and  are  solved  at 
each  time  step  to  determine  the  potentials  from  the  charge  and  current 
densities.  The  electric  and  magnetic  fields  are  then  obtained  from  the 
potentials  in  the  usual  way.  The  field  equations  are  differenced 
implicitly  in  both  the  radial  and  axial  directions,  but  an  alternate 
version  with  explicit  differencing  in  the  axial  direction  has  also  been 
used.  Once  the  fields  are  calculated,  the  code  pushes  all  particles 
relativistically  in  Cartesian  coordinates  using  standard  techniques.  In 
order  to  determine  the  densities,  we  use  a  particle  lay-down  scheme  which 
is  quadratic  in  the  radial  direction  and  linear  (or  Nearest  Grid  Point  in 
the  alternate  version)  in  the  axial  direction.  The  code  should  be  valid 
for  moderate  energy  as  well  as  high  energy  beam  propagation  studies. 
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IX 


Table  1 


Erosion  rates  from  simulation  results  and 
predicted  by  Eq.  (7). 


<kA) 

f 

P0-e  ( FRIEZR) 

2 

0.2 

0.024 

2 

0.5 

0.039 

2 

1.0 

0.094 

5 

0.2 

0.043 

5 

0.5 

0.085 

5 

1.0 

0.16 

10 

0.2 

0.061 

10 

0.5 

0.15 

10 

1.0 

0.28 

corresponding  values 


eo-ep  [Eq.  (7)] 

0.022 

0.046 

0.084 

0.040 

0.090 

0.17 

0.062 

0.14 

0.26 
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Fig.  1  —  E  vs  f.  R]r  vs  f,  and  nt  and  nh  vs  f  arc  plotted  separately  at  <  t  =  7  meters  tor  a  beam  with 
parameters:  !h  =  5  kA.  /  =  1.0.  X„  =  10.  rh  =  1  cm.  r,  =  2  cm  and  r„  =  33  cm 
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Fie.  4  /.  vs  R ,  ,  ■ 

the  same  as  in  Fiji.  3. 

vs  j-,  and  n,  and  nh 

vs  f  arc  plotted  separately  at  ct  =  30  meters. 

The  parameters  are 

cm 


8x10 
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6x1010 
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Fig.  4  (Continued)  E.  vs  f,  R r2  vs  and  n,  and  nh  vs  f  are  plotted  separately  at  ct  =  30  meters.  The 
parameters  are  the  same  as  in  Fig.  3. 
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Fig.  5  —  /cn, |  vs  f  plotted  at  cr  =  4  and  40  meters  shows  an  increase  in  the  rise  length,  fr.  and  the 
development  of  a  low-emittance  population  at  the  beam  head.  Areas  with  low  transverse  temperature, 
7  <  T„/ 2.  are  shaded.  Parameters  are  lh  =  1  kA,  /  =  1,  X„  =5  and  rh  =  r,  =  2  cm.  E.  =  0  is 
specified  in  the  particle  push  routine. 


(cm) 


Fig.  5  —  (Continued)  /cnc|  vs  f  plotted  at  c t  =  4  and  40  meters  shows  an  increase  in  the  rise  length,  fr.  and 
the  development  of  a  low-emittance  population  at  the  beam  head.  Areas  with  low  transverse  temperature, 
T  <  T„/l,  are  shaded.  Parameters  are  Ih  =  I  kA,  /  =  1,  X„  =5  and  rh  =  r,  =  2  cm.  E.  =  0  is 
specified  in  the  particle  push  routine. 
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Fig.  6  —  The  coordinate  location  of  the  T  =  T„f2  point,  fr,  is  plotted  vs  rr  for  the  simulation  pictured  in 
Fig.  5.  As  the  asymptotic  near-equilibrium  state  is  reached,  this  point  will  become  stationary  with  respect  to 
the  beam  body,  attaining  a  slope  d{r/dz  =  d£h/dz  =  1-/3,  where  d{h/dz  is  the  velocity  of  the  beam  bod) 
relative  to  the  simulation  coordinate  frame  and  z  -  ct.  A  line  with  slope  1  -  0  is  also  plotted. 
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